Mechanical Control of Spin States in Spin-1 Molecules 
and the Underscreened Kondo Effect* 
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Controlled stretching of individual transition-metal complexes enables direct manipulation of the 
molecule's spin states and tests of predictions for the underscreened S = 1 Kondo effect. 

The ability to make electrical contact to single molecules creates opportunities to examine fun- 
damental processes governing electron flow on the smallest possible length scales. We report ex- 
periments in which we controllably stretch individual cobalt complexes having spin S = 1, while 
simultaneously measuring current flow through the molecule. The molecule's spin states and mag- 
netic anisotropy were manipulated in the absence of a magnetic field by modification of the molecular 
symmetry. This control enabled quantitative studies of the underscreened Kondo effect, in which 
conduction electrons only partially compensate the molecular spin. Our findings demonstrate a 
mechanism of spin control in single-molecule devices and establish that they can serve as model 
systems for making precision tests of correlated-electron theories. 



The electronic states of atoms and molecules depend 
on the symmetry imposed by their local environment. A 
simple manifestation occurs in the attachment of ligands 
to a transition-metal ion to form a coordination com- 
plex, which breaks spherical symmetry and causes split- 
tings within the ion's initially degenerate d-orbitals. For 
a complex having spin 5 > 1, additional distortions of 
the ligands combined with spin-orbit coupling can cause 
splittings within the initially (25 + l)-degenerate spin 
states, giving rise to magnetic anisotropy [1, 2]. To 
study the effects of symmetry-breaking distortions, we 
stretched individual 5 — 1 molecules within mechanically 
controllable devices [3]. Simultaneous electron transport 
measurements showed that stretching lifts the degener- 
acy of the 5 = 1 ground state and enables control of the 
magnetic anisotropy. The same devices also enabled tests 
of predictions [4] for the temperature dependence of the 
underscreened 5 = 1 Kondo effect [5-7]. 

We studied the Co(tpy-SH) 2 complex (Fig. 1A), in 
which a Co ion resides in an environment of approxi- 
mately octahedral symmetry, through its coordination to 
six N atoms on two terpyridine ligands. When attached 
to gold electrodes and cooled to low temperature, this 



molecule exhibits the Kondo effect [8]; the spin of the 
molecule is screened by electron spins in the electrodes, 
leading to a peak in the electrical conductance at zero 
bias voltage, V [9]. We used the Kondo effect as a spec- 
troscopic probe to interrogate the molecular spin [10]. 
Previous electron-transport studies of individual metal 
complexes have probed rich behavior including not only 
Kondo physics [8, 11-13], but also vibrational excitations 
[13-15] and molecular magnetism [16, 17]. 

We used mechanically controllable break-junction de- 
vices (Fig. IB) to stretch individual molecules while mea- 
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FIG. 1: Single-molecule electrical devices based on a transition- 
metal complex. (A) Chemical structure of Co(tpy-SH)2, where 
tpy-SH = 4'-mercapto-2,2':6',2"-Terpyridine. (B) Schematic of the 
mechanically controllable break junction. (C) Scanning electron 
micrograph of a suspended break junction in false color, prior to 
molecule deposition. 
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FIG. 2: Splitting of the Kondo peaks as a function of mechanical 
stretching. (A),(C) dl /dV as a function of bias voltage and elec- 
trode spacing at T = 1.6 K for (A) device A and (C) device B. 
(B),(D) Line cuts showing dl/dV as a function of V for different 
values of molecular stretching for (B) device A and (D) device B. 
(E) Temperature dependence of the conductance for device A on 
the stretched side of the splitting transition. (F) Breaking of the 
degeneracy among the S = 1 triplet states by uniaxial distortion. 



suring their conductance. The devices were made by fab- 
ricating a Au wire suspended above a thin Si substrate 
(Fig. 1C), depositing molecules, and then using electro- 
migration [18] to create a molecular-scale break in the 
wires before beginning studies in which we tuned the 
distance between electrodes by bending the substrate. 
Details, statistics, and discussion of control experiments 
are provided in the Supporting Online Material, section 
SI [19]. We have also measured the same Co(tpy-SH) 2 
complex using fixed Au electrodes [8, 11, 13, 15-17]. 

Measurements of differential conductance (dl/dV) as a 
function of increasing electrode spacing at a temperature 
T = 1.6 K are shown in Fig. 2. At the initial position of 
the electrodes, the dl/dV spectra exhibited a single peak 
centered at V — with amplitude of order (but less than) 
2e 2 /h, a signature of Kondo-assisted tunneling through 
the molecule. As we stretched each molecule, the single 
conductance peak split into two beyond a value for the 
change in electrode spacing that varied from device to 
device. For the two devices displayed in Fig. 2 and for two 
others (see section S2 [19]), we could reproducibly cross 



this transition back and forth between one peak and two. 
For the stretched molecules (Fig. 2E), the temperature 
dependence of dl/dV at V = showed a nonmonotonic 
dependence similar to the ^-dependence upon increasing 
|V| from 0. 

The observed splitting of the Kondo peak as a func- 
tion of stretching is in striking contrast to a previous 
study of the spin-| Kondo effect in Ceo molecules [20], in 
which varying the electrode spacing modified the height 
and width of the Kondo resonance but did not split the 
peak. We show that the peak splitting for the Co(tpy- 
SH) 2 complex is caused by a higher-spin 5 = 1 Kondo 
effect, together with the breaking of degeneracy within 
the 5 = 1 triplet ground state caused by molecular dis- 
tortion [2] . For an unstretched 5=1 ion in a ligand field 
with octahedral symmetry, the triplet states are strictly 
degenerate according to group theory. However, if the 
molecule is stretched axially (the z-axis), the S z = 
state will be lowered by a zero-field splitting energy D 
below the 5 2 = ±1 states, corresponding to a uniax- 
ial spin anisotropy (Fig. 2F). This broken degeneracy 
quenches the Kondo resonance near V = and causes 
conductance peaks at V — ±D/e due to inelastic tunnel- 
ing. The distortion-induced breaking of the ground-state 
degeneracy for an 5 = 1 molecule is in contrast to the 
physics of half-integer spins [21] for which time-reversal 
symmetry mandates that the ground state remains a de- 
generate Kramers doublet. To establish this picture of 
stretching-induced control of spin states, we show that 
the spin of our molecules is indeed 5=1, based on tem- 
perature and magnetic-field studies, and we verify that 
stretching produces spin anisotropy by measuring how 
the level splitting depends on the direction of an applied 
magnetic field. 

The temperature dependence of the Kondo conduc- 
tance is predicted to depend on the molecular spin 5 and 
the number of screening channels in the electrodes. To be 
fully screened at zero temperature, a molecule with spin 
5 requires coupling to 25 screening channels [5]. In our 
experimental geometry, there are two screening channels 
consisting of linear combinations of states from the two 
electrodes with different couplings to the impurity J\ and 
Ji that result in two Kondo temperature scales Tki and 
Tki [22] (we assume Tki < Tk%)- The Kondo temper- 
atures depend exponentially on the couplings, so in the 
typical case Tki « Tk2- If 5 > \, the result is that 
over the range Tki << T < Tk2 the original spin is only 
partially screened, to a value 5— \. Relative to the fully 
screened 5 = \ Kondo effect, this underscreened effect 
should produce a much slower rise in the conductance as 
T decreases below Tki [4, 6, 7]. Henceforth we denote 
Tki as simply Tk ■ 

Figures 3A and 3B show the measured zero-bias con- 
ductance G(T) for the unstretched configuration of de- 
vices A and B with fits to numerical renormalization 
group (NRG) predictions for the fully screened 5 = \ 
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FIG. 3: Temperature dependence corresponding to the under- 
screened Kondo effect. (A),(B) Fits of zero-bias conductance as 
a function of temperature for (A) device A and (B) device B to 
NRG calculations for S = i, 1, and |. (C) Separate fits to the 
S = „-, 1] an d I Kondo models, combining data from seven de- 
vices. The collapse of the data on the predicted scaling curve for 
the underscreened 5 = 1 Kondo model is excellent; the normalized 
root-mean-square scatter for the S = | ansatz is greater by 44%, 
and for S = § by 320%. For the S = 1 fits, the seven Kondo 
temperatures are 1.1 K, 28 K, 62 K, 81 K, 84 K, 180 K, and 210 K. 



Kondo model [23, 24] and the underscreened S — 1 and 
| models [4] (for fitting details, see section S3 [19]). 
Each fit has two adjustable parameters, Tk and the zero- 
temperature conductance G(0). We find that G{T) devi- 
ates strongly from the form for the S = \ Kondo effect, 
and instead agrees quantitatively with the prediction for 
the S = 1 underscreened case. We performed the same 
analysis for 10 unstretched Co(tpy-SH)2 devices (section 
S4 [19]). Seven showed similarly unambiguous agreement 
with underscreened Kondo scaling, and all of these gave 
superior fits to the S — 1 prediction than to S = |. Dis- 
cussion of the other 3 samples is provided in section S4. 

In Figure 3C, for the 7 devices with underscreened 
characteristics, we plot G(T)/G(0) versus T/Tk using 
the parameters G(0) and Tk obtained from separate fits 
to the S = |, 1, and | models. This allows us to test how 
well the data can all be described by the predicted scaling 
curves. We used the root mean square deviation of the 
data from each theory curve normalized to the average 
of the scaled data as a goodness-of-fit metric. The S — 
1 fit is best: the deviation of points for the S = J ansatz 
is greater by 320%, for S = | by 44%, and for larger 



S by more than 80%. Additional confirmation that S 
= 1 is discussed below based on magnetic field studies. 
An independent [25] observation of the temperature scal- 
ing predicted for the underscreened Kondo effect has also 
been reported recently for a single S = 1 Cgo device [26] . 

Our interpretation that the stretching-induced Kondo 
splitting is caused by the breaking of spin degeneracy 
is confirmed by the presence of spin anisotropy in the 
stretched state - the magnetic field (B) dependence of 
the Kondo splitting differs depending on the angle at 
which B is applied relative to the stretching axis. Fig- 
ure 4A shows the expected B-dependences of the energy 
levels in an S = 1 multiplet for D > 0, for several field an- 
gles, and Fig. 4B displays the energy differences between 
the first excited state and the ground state, which deter- 
mine the Kondo peak positions. For a field orientation 
approximately perpendicular to the anisotropy axis (9 = 
90°), when gfisB « D the i?-dependence of the peak 
position should be very weak, while when gfisB grows 
large enough to be comparable to D, the peak should 
shift gradually to larger values of |V| with positive cur- 
vature. For B parallel to the anisotropy axis (6 = 0°), 
the -B-dependence of the peak position should be linear 
and much stronger. Our measurements show a depen- 
dence on the orientation of B just as expected within 
this model. Figures 4C and D show the i?-dcpendencc 
for different degrees of stretching for device A, with 9 « 
90°. At a relatively large degree of stretching (Fig. 4C), 
D = 3.45 meV was larger than g^sB at our largest field 
(~1.5 meV), and the Kondo peak positions were almost 
independent of B. The slightly negative slope d\V\/dB 
(corresponding to an effective g-factor of 0.21 ± 0.01) 
suggests that the stretching axis is not exactly perpen- 
dicular to B. For a smaller degree of stretching (Fig. 
4D), gfisB was comparable to D at large fields, and 
the Kondo peaks shifted weakly to larger |V| as a func- 
tion of B with positive curvature. Figures 4E and F 
show data for B approximately parallel to the molecu- 
lar axis for fixed-electrode devices, which exhibit zero- 
field splitting without intentional mechanical stretching 
(among devices that had Kondo features, ~10% of the 
fixed-electrode devices and ~20% of the adjustable de- 
vices exhibited a split Kondo peak at the initial electrode 
spacing, suggesting that these molecules happened to be 
strained initially). In Figures 4E and F, the Kondo peaks 
shift linearly with B, with much larger slopes d\V\/dB 
than for 9 sa 90°, corresponding to ^-factors of 2.6-2.7 (± 
0.3). In device D, the peaks shift initially to smaller |V|, 
pass through zero, and then move to larger \V\ (Fig. 4F). 
From the negative slopes of d\V\/dB in Figures 4C and 
4F, we can confirm that D > 0. The scale of splittings we 
measure, ~0-5 meV, agrees with our simulations (section 
S6 [19]) and is typical of zero-field splittings in distorted 
coordination complexes [2]. 

These magnetic field data provide further evidence 
that the spin state of the molecule is S = 1, rather than 
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FIG. 4: The magnetic field evolution of the Kondo peaks demonstrates the presence of spin anisotropy. (A) Energy eigenvalues (for 
S = 1) of the model spin-anisotropy Hamiltonian H = — 9A*B ■ S + DS^ for 4 different field angles with respect to the anisotropy axis. 
The red arrows indicate the lowest-energy inelastic transitions corresponding to the finite-bias Kondo peaks. (B) Predicted Kondo peak 
position versus magnetic field at several field angles, with energies normalized to the zero-field splitting. Inset: Kondo peak splitting (for 
g = 2) as a function of magnetic field for all angles in the absence of spin anisotropy. (C),(D) Kondo peak positions versus magnetic 
field at T = 1.6 K in device A for two different values of stretching, with the magnetic field oriented approximately perpendicular to the 
molecule axis. The data in (D) can be fit for any field angle between 70° and 90°, giving parameters in the range D = 2.0 - 2.2 meV and 
g = 2.5 - 2.8. (E),(F) The magnetic field dependence in fixed-electrode devices C and D with the magnetic field approximately along the 
molecular axis. The insets contain color plots of dl/dV versus bias voltage and magnetic field. 



|, |, or any half-integer. For half-integer spins with D > 
0, the ground state for B — is the degenerate Kramers 
doublet S z — ±i. Regardless of the direction of B or 
the degree of stretching, this doublet will give rise to a 
Kondo peak that shifts approximately linearly with B 
and extrapolates to zero splitting at B — [21], in con- 
tradiction to the data (section S7 [19]). Integer values of 
S = 2 can be ruled out based on the temperature scaling 
described above and because previous studies of six-fold 
coordinated Co complexes have found almost exclusively 
S < |, with exceptions only for fluoride ligands [1, 2]. 
We conclude that only S = 1 can explain the measure- 



ments, and from this we identify the charge state of the 
metal center to be Co 1+ (section S8 [19]). 

In coordination chemistry, the existence of zero- 
field splittings induced by molecular distortion is well- 
established, but the ability we demonstrate to contin- 
uously distort an individual molecule while simultane- 
ously measuring its zero-field splitting opens the possi- 
bility for dramatically more detailed and precise compar- 
isons with theory. For correlated-electron physics, our re- 
sults demonstrate that single-molecule electrical devices 
can provide well-controlled model systems for studying 
5 = 1 underscreened Kondo effects not previously realiz- 



able in experiment. Our work further demonstrates that 
mechanical control can be a realistic strategy for manip- 
ulating molecular spin states, to supplement or replace 
the use of magnetic fields in proposed applications such 
as quantum manipulation or information storage [27, 28]. 
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To make the mechanically controllable break-junction devices, we start by fabricating con- 
tinuous gold lines (32 nm thick, 500 nm long, and with a 50 nm-wide constriction) suspended 
40 nm above a 200 fim-thick Si wafer [S1,S2]. Fabrication of the fixed-electrode devices follows 
procedures described in Refs. [S3,S4]. To incorporate the molecules into either type of device, we 
clean unbroken wires in an oxygen plasma to remove organic contaminants and then immerse the 
samples into a <0.1 mM solution of Co(tpy-SH)2 with PF~ counterions in acetonitrile, allowing 
the thiol end groups of the molecules to attach to the gold. We synthesized the molecules by 
the process described in Ref . [S3] and purified them by flash chromatography using an alumina 
column, yielding a final purity >99% prior to deposition. After the sample chip is removed from 
the molecular solution, excess solution is blown off with nitrogen gas and the samples are cooled 
to low temperature (T = 1.6 K). We then use electromigration [S5] to create a molecular-scale break 
in the wires. The electrode motion is calibrated [S6] from the tunneling conductance of junctions 
lacking any added molecules. 

After electromigration, we observed a peak in dl/dV at V = 0, characteristic of the Kondo effect, 
in approximately 14% of ~250 mechanically controllable devices and 18% of ~200 fixed-electrode 

1 



devices made with Co(tpy-SH)2 solution. Approximately 3% of devices exhibited a high-resistance 
Coulomb blockade characteristic without a Kondo feature, and the rest had featureless dl/dV vs. 
V traces indicative of simple tunneling. Among the devices with a conductance feature near zero- 
bias, at the initial electrode spacing ~20% of the adjustable devices and ~10% of the fixed-electrode 
devices had split peaks (suggesting that the molecule is already stretched); the peak positions 
ranged up to \V\ = 5 mV. Four adjustable devices possessed sufficient stability to repeatedly 
observe the transition between split and unsplit peaks as a function of stretching. No devices 
showed a decrease in the Kondo splitting with stretching; two showed a zero-bias Kondo peak 
that did not split with stretching (see Section S5). Among 47 adjustable and 24 fixed-electrode 
devices prepared as control samples with acetonitrile alone, none contained a zero-bias peak in 
dl/dV and only one exhibited Coulomb blockade. 



S2. Additional devices for which the Kondo peak splits with stretching 

In Fig. SI we show data for two devices (Devices I and J), in addition to the two described 
in the main text, that display a splitting of the Kondo peak with stretching. While we were able 
to tune from a single Kondo peak to split Kondo peaks and then return to a single peak on a 
subsequent reverse scan with Devices I and J, these junctions were not sufficiently stable to make 
detailed temperature and magnetic field measurements. 

B 
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Figure SI: Two additional devices that showed a splitting of the Kondo peak with stretching. 



S3. Fitting to Numerical Renormalization Group (NRG) predictions for G(T) 



NRG calculations for Kondo-assisted tunneling via an impurity of spin S coupled to a single 
screening channel predict the temperature dependence of the conductance in the form of a discrete 
set of points, G/G(0) as a function of T [S7]. The calculations were carried out with a logarithmic 
discretization parameter A = 1.5 [S8] and retaining 1700 states per NRG iteration. In order to 
fit experimental data to these NRG predictions, we first approximate the NRG results by a set 
of analytical fitting functions having the scaling form Gg(T) = G(0) /s(T/Tk), where Tk is defined 



so that G(Tk) = G(0)/2, as was done in Refs. [S8, S9]. We do not wish to imply any physical 
significance to the form of these fitting functions; they are merely an attempt to approximate the 
NRG results with a minimum number of adjustable parameters. We use as a starting point the 
phenomenological Goldhaber-Gordon form [S9] for the temperature dependence of the spin-1/2 
Kondo effect 



gs(T/T K ) = 



11^ 

1 K 



-as 



(SI) 



but we relax the condition Es - 2 (applicable for a fully screened Kondo effect) since for S > 1/2, 
the ground state of an underscreened Kondo system is not a Fermi liquid. Note that the quantity 
T' K of Eq. (SI) is related to the actual Kondo temperature by 



T = Tk 

k (2 1 /"s - l)V£s ' 



(S2) 



We find that as the temperature is reduced, functions of the form in Eq. (SI) saturate more quickly 
to the maximum value of 1 than do the NRG results for the underscreened Kondo effects - the 
calculated conductance is not yet saturated at the lowest temperatures T/Tk ~ 10~ 4 available in 
the results. To better account for this slow saturation, we find over many decades of temperature 
around Tk that good fits to the NRG results can be obtained by including a small offset parameter 
As, such that for our phenomenological scaling functions we use the form 



GUT) = G(0) [(1 - A s )gs(T/T K ) + A s /2] . 



fs(T/T K ) 



(S3) 
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Figure S2: Fits of NRG conductances G(T) to Eq. (S3) for S = 1/2, 1, 3/2, 2, and 5/2. 



In the above, the coefficient (1 - As) helps to account for the incomplete saturation of conductance 
for underscreened Kondo effects at the lowest temperatures available in the NRG results, and the 
second term Ag/2 is added to retain consistency with the definition of the Kondo temperature; one 
can verify by inspection that the term allows the condition G* s (Tk) - G(0)/2 to be satisfied. The 
reader should not be confused by the fact that G* s (0) + G(0), because in our fits we will be concerned 
only about temperatures within several decades of Tk, not with the ultimate low-temperature limit. 

In Fig. S2, we show fits of Eq. (S3) to the NRG data for S - 1/2, 1, 3/2, 2, and 5/2, where 
£,$/ oc$, arid Ag are used as fitting parameters. The case of S = 1/2 is an exception - here we fixed 
£s = 2 and A$ — as was done in Ref. [S9]. The fitting parameters are summarized in Table SI. 
The NRG results are very well described by Eq. (S3) for over 7 decades in T/Tk, with particularly 
good agreement in the regime 10~ 3 Tk < T < Tk, which is the most relevant for comparison to 
our experiments. In order to fit the predictions of the Kondo models for different spins S to our 
experimental data, we employ Eq. (S3) keeping fixed the values of the parameters £,s, &s, and A$ 
as listed in Table SI, and use G(0) and Tk as the only adjustable fitting parameters. In fitting to the 
experimental data, we did not include a background term as an adjustable parameter in Eq. (S3). 



Table SI: Fitted parameters determined by using Eq. (S3) to approximate the NRG predictions for 
the temperature dependence of conductance for a Kondo impurity of spin S with a single 
screening channel (giving a fully screened Kondo model for S — 1/2 and underscreened 
forS> 1/2). 



Spin 


£s 


a s 


A s /2 


S = l/2 


2 


0.220 + 0.005 





S = l 


0.745 + 0.009 


0.506 + 0.009 


0.041 + 0.001 


S = 3/2 


0.483 + 0.004 


0.670 + 0.008 


0.062 + 0.002 


S = 2 


0.349 + 0.004 


0.912 + 0.013 


0.093 + 0.002 


S = 5/2 


0.288 + 0.004 


1.116 + 0.023 


0.124 + 0.003 



S4. Fitting the measured G(T) to Kondo scaling predictions for different values of S, 
for all samples 

We measured the detailed temperature dependence of the unsplit Kondo signal in ten Co(tpy- 
SH)2 devices, including two devices with adjustable electrodes (A and B) and eight with fixed 
electrodes. For each set of data GCT) we performed fits to the predicted scaling curves (Eq. (S3)) 
for the fully screened S - 1/2 Kondo model and the underscreened S — 1 and 3/2 Kondo models, 
with two fitting parameters for each fit, the zero-temperature conductance G(0) and the Kondo 
temperature Tk- (To be clear, different best-fit values of G(0) and Tk were determined for each 
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B. S= 1/2 scaling for two Co(tpy-SH) devices (likely a different charge state of the molecule) C. Ambiguous scaling 
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D. S= 1/2 scaling from previously studied C 60 devices 
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Figure S3: Fits of G(T) to Eq. (S3) for ten Co(tpy-SH)2 devices. (A) Seven devices show a tem- 
perature scaling that deviates from the S - 1/2 fully screened Kondo prediction, and is 
instead in good agreement with the underscreened S - 1 Kondo effect. (B) Data from 2 
devices are consistent with a spin-1/2 Kondo effect, likely arising from a molecule with 
a Co 2+ charge state. (C) One device did not unambiguously fit any of the predicted 
forms for a single Kondo impurity. (D) S - 1/2 scaling observed in previously studied 
C6o devices [S2]. 



value of S.) The results are shown in Fig. S3 for all ten unstretched Co(tpy-SH)2 devices, seven of 
which showed S — 1 scaling, two of which showed S - 1/2 scaling, and one device that showed 
ambiguous scaling. For comparison, we also show the same analysis for three previously studied 
Cgo devices with S - 1/2 [S2]. 

The first conclusion we can draw is that generally there is no ambiguity in distinguishing 
between realizations of the fully screened S - 1/2 Kondo effect and the underscreened Kondo effects, 
because the underscreened models predict a very slow approach to saturation with decreasing 
temperature, whereas the fully screened models saturate rapidly at low T. All of the devices we 
have measured, except for the one exception noted, give excellent fits to either the fully screened 
S - 1/2 curve or the underscreened curves, with a very poor fit to the other. From this we conclude, 
first, that the seven devices shown in Fig. S3A have S > 1. That the temperature dependence for 
each of these samples can be fit well with a single value of Tk indicates that we measure current 
flow via a single molecule in each device. 

Distinguishing between the different models of the underscreened Kondo effect (e.g., S = 1 vs. 
3/2) based solely on the temperature dependence is more challenging because the scaling curves 
are similar, with most of the difference due to a few points at the high and low temperature ends of 
the fitting ranges. However, all seven of the samples that exhibit an underscreened characteristic 
exhibit better agreement with the S — 1 curve than with S - 3/2 or higher spin. As discussed 
in the main text, the magnetic field dependence of the Kondo peak splitting provides additional 
confirmation that S - 1. 

For the devices in Fig. S3B that show S — 1/2 scaling, our interpretation is that the molecule 
in these devices likely retained a different charge state (e.g., Co 2+ ) than the molecules exhibiting 
the underscreened Kondo effect. The presence or absence of nearby counterions or other charged 
impurities, or differences in local work functions for the disordered gold electrodes, could easily 
lead to large differences in the local electrostatic potential at the position of the molecule through 
which current flows, thereby providing in effect a local gate to change its charge state. The Co 2+ 
state should have a ground state spin S = 1/2 at low temperature [S10], which can produce the 
conventional fully screened S — 1/2 Kondo effect. 

The one device showing a temperature dependence different from both the S - 1 and S - 1/2 
predictions (Fig. S3C) may possibly be explained by conduction through two molecules in parallel, 
with different Kondo temperatures or possibly to being in a mixed-valence, rather than a purely 
Kondo, regime [Sll]. 



S5. Two Co(tpy-SH) 2 devices did not show a splitting of the Kondo peak with stretching 

Increasing the electrode spacing did not result in a splitting of the Kondo peak in two devices 
that we measured. The dl/dV of one such device as a function of increasing electrode spacing is 
shown in Fig. S4. Stretching modifies the height and width of the Kondo peak, but the peak does 
not split. This behavior might be explained in several ways. One is that the molecule may not be 
strongly bonded to both electrodes so that electrode displacement does not significantly stretch the 
metal complex. The effect of electrode displacement in this case would be only to reduce the Kondo 
conductance and temperature without significantly affecting the symmetry of the molecule [S2] . 
A second possibility is that the spin state of the molecule for these two devices is S = 1/2 (i.e., in 
the Co 2+ state), as discussed above. 
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Figure S4: dl/dV traces with increasing electrode spacing for a device at T = 1.6 K, which did not 
show a splitting of the Kondo peak with stretching. 



S6. Calculations of zero-field splittings 



We carried out calculations of the zero-field splitting for the isolated molecule spin-triplet 
state within the DFT coupled-perturbed spin-orbit coupling formalism [S13], as implemented in 
ORCA [S14]. Precise reproduction of the observed experimental D values is challenging because 
of the dependence on the molecular environment. It is sensitive in particular to the overall spin 
density on the cobalt, which gives some insight into the mechanism of how the variation in D may 
be affected at a microscopic level. For the molecule at equilibrium geometry, we find that the zero- 
field splitting is ~0.1 meV, consistent with the expectation of a small value for nearly octahedral 
symmetry. As we tune the exchange contribution in the density functional to shift the spin density 
onto the cobalt atom, the zero-field splitting can increase. For an exchange contribution that 
produces a similar magnetic moment (1.2 Bohr magnetons) on the cobalt atom as observed in LDA 
+ U calculations on the molecule-electrode system, we obtain a D of 2.7 meV (using 80% exchange). 



This is of the same sign and similar in magnitude to the splittings observed in the experiment. The 
amount of exchange required is quite high, but the primary effect (to shift spin density onto the 
cobalt atom) might arise in the experiment from attachment to the electrodes. 

We also tried to obtain the dependence of D on molecular geometry. This was complicated 
by numerical artifacts arising from the solution of the coupled-perturbed equations treatment. 
Redistribution of the unpaired electrons between the Co ion and the ligands might also contribute 
to changes in the magnitude of D as the metal-ligand coupling changes. However, variations on 
the order of 1 meV were calculated. 



S7. Spin energy levels as a function of magnetic field for S = 3/2 



In Figures 4A and 4B of the main manuscript, we plotted the spin energy levels and the 
Kondo peak position as a function of magnetic field for several different field angles relative to the 
stretching axis z, for a simple spin Hamiltonian H = -g/igB • S + DSl with S — 1 in the presence 
of a spin anisotropy energy D > 0. In Fig. S5 we plot the results of the same model Hamiltonian 
for S — 3/2. For any half-integer spin value with D > 0, the ground state for B — is a degenerate 
Kramers doublet with S z — +1/2. When a magnetic field is applied, no matter what the orientation 
of B or the magnitude of D, the Kramers doublet splits approximately linearly (for gfiBB < D), and 
therefore the Kondo peak position would always undergo an increase that extrapolates to zero for 
B - 0. This is in contradiction to the field dependence that we measure for stretched Co(tpy-SH)2 
devices, allowing us to rule out the possibility that S = 3/2, or any other half-integer value. 
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Figure S5: Magnetic field evolution of S - 3/2 energy levels in the presence of anisotropy. (A) 
Energy eigenvalues of the model spin-anisotropy Hamiltonian H = -gjUgB • S + DS^ for 
4 different field angles with respect to the anisotropy axis. The red arrows indicate the 
lowest-energy inelastic transitions corresponding to the finite-bias Kondo peaks. (B) 
Predicted Kondo peak position versus magnetic field at several field angles. 



S8. Identification of the charge state of the measured Co(tpy-SH) 2 

We can identify the charge state of the measured Co(tpy-SH)2 complex based on the deter- 
mination that the ground-state spin is S = 1. For a complex in equilibrium with counterions, the 
metal center would be in the Co 2+ state with an electronic configuration of 3d 7 and spin S — 1/2 or 
3/2 with predominantly S - 1/2 at low temperature for this complex [S10], inconsistent with the 
measured ground state of S = 1. Instead, charge transfer of an electron from the gold electrodes to 
the complex must result in a 3d 8 Co 1+ state, which has an S = 1 ground state for an approximately 
octahedral ligand field (Fig. S6A). The other possible low-energy even-electron charge state, Co 3+ , 
has S - for 6-fold coordinated complexes [S12]. 

free Co 1+ ion octahedral 

(3 c/ 8 configuration) field 

t t 



Mltlt 

3dorbitals 
S=1 



ft&tt 
L 



Figure S6: Electronic configuration of the cobalt complex. For a Co 1+ ion coordinated to ligands 
with approximately octahedral symmetry, the 5 metal d-orbitals split into a set of tig 
and eg orbitals, and a ground state with S - 1 is expected. 



S9. Consideration of alternate mechanisms for the stretching-induced Kondo peak 
splitting 

We have considered and can rule out two potential alternative mechanisms for the splitting 
of the Kondo peak induced by stretching. The "two-stage" S = 1 Kondo effect [S15,S16], relevant 
for temperatures less than the Kondo energy scales for both screening channels (Tki and Tjq), 
can lead to a non-monotonic bias and temperature dependence similar to what we observe in 
the stretched molecules. However, for this mechanism to explain our data both Tki (reflected in 
the overall width of the Kondo peak) and Tki would have to increase with stretching. This is 
unphysical because stretching should weaken the coupling of the metal center to the electrodes. 
A splitting of an S = 1 Kondo signal could also result from a Jahn-Teller mechanism if distortion 
caused a non-degenerate S = singlet to shift below the S = 1 triplet to become the ground state 
of the molecule. A triplet-singlet crossing has been suggested as the explanation for gate-voltage- 
induced splittings of Kondo peaks associated with an even number of electrons in GaAs [SI 7], 
carbon nanotube [S18], and molecular [S19, S20] quantum dots. However, ab initio calculations 



indicate that for the range of stretching in our experiment, the S - 1 triplet state is always the 
ground state of Co(tpy-SH)2 (see Section S10). Furthermore, neither a two-stage Kondo effect nor 
a triplet-singlet crossing would explain the magnetic-field anisotropies that we measure. 



S10. Calculations of the triplet-singlet energy gap 

We investigated the singlet and triplet energetics of the Co(tpy-SH)3, + complex with and 
without attachment to electrodes. Isolated molecule calculations were performed using the ORCA 
and MOLPRO packages [S14,S21]; molecule-electrode calculations were performed using the VASP 
package [S22]. We simulated stretching of the molecule in two ways (i) altering the S-S distance, 
allowing geometric relaxation of the other atoms, and (ii) moving the axial pyridyl moieties in a 
rigid fashion (see Fig. S7). 

Using density functional theory (B3LYP/Def2-TZ2P), we computed the closed-shell (spin- 
unpolarized) triplet-singlet gap. This tests the Jahn-Teller type scenario mentioned in the previous 
section - whether mechanical distortion might cause an orbital energy splitting favoring the singlet 
over the triplet state. We found, however, that both types of stretching lead to the triplet state 
being below the singlet state by about 0.4 eV (Fig. S8), much larger than any energy scale ob- 
served in the experiment. We also carried out multireference complete active space self -consistent 
field (CASSCF/6-31G) calculations using a 14-electron 14-orbital active space and this yielded a 
comparable singlet-triplet gap (0.55 eV) near the equilibrium geometry. 

Electrodes: To test the effect of attaching to electrodes, we carried out periodic density func- 
tional calculations (PAW, 350 eV plane-wave cutoff) with the molecule inserted between two model 
gold electrodes. Here, pure functionals such as LDA and PBE did not yield magnetization of the 
molecule in the triplet state, an artifact of the density functional approximations. Moving to the 
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Figure S7: Two simulations of molecule stretching. (A) Stretching by increasing the sulfur-sulfur 
distance while relaxing the positions of the other atoms, and (B) stretching by rigidly 
displacing the axial pyridyls. 
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Figure S8: Calculated triplet-singlet energies as a function of stretching obtained by (A) varying the 
sulfur-sulfur distance (with relaxation of the other atoms), and (B) rigid displacement 
of the pyridyl groups. Energies are measured relative to the lowest triplet energy. 



LDA + U formalism (Co: U = 8.0 eV, / = 1.0 eV) [S23] recovered a magnetic cobalt atom in the triplet 
state, with a corresponding singlet-triplet gap of 0.52 eV at a near-equilibrium molecule-electrode 
geometry. This appears consistent with our isolated molecule studies. 

Other low-lying states: We also investigated the presence of other low-lying states in the 
spectrum. One possibility is an antiferromagnetic (open-shell) singlet where the cobalt atom is 
formally Co(III) and the two terpyridine ligands are formally negatively charged with an unpaired 
electron, antiferromagnetically coupled to each other. In the isolated molecule case, we found the 
AFM singlet to be lower in energy than the closed-shell singlet. The AFM singlet-triplet gap is 
very sensitive to the precise theoretical treatment. By varying the amount of exact exchange in the 
density functional, 0% (BP86), 20% (B3LYP), 50% (BHLYP) we observed gaps of 19.8 meV, 3.1 meV, 
1.9 meV. In the case of attached electrodes, we did not find an AFM singlet solution using the LDA 
and PBE functionals, likely due to the same self-interaction errors in the density functionals that 
led to the disappearance of the magnetism in the triplet state mentioned above. Nonetheless, we 
consider it likely that the AFM singlet state may be a low-lying state in these systems. 
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